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One way to select a preferred frame from gravitational radiation is via the principal axes of 
an average of the action of rotation group generators on the Weyl tensor at asymptotic 
infinity. In this paper we evaluate this time-domain average for a quasicircular binary using ap- 
proximate (post-Newtonian) waveforms. For nonprecessing unequal-mass binaries, we show the 
dominant eigenvector of this tensor lies along the orbital angular momentum. For precessing bina- 
ries, this frame is not generally aligned with either the orbital or total angular momentum, working 
to leading order in the spins. The difference between these two quantities grows with time, as the 
binary approaches the end of the inspiral and both precession and higher harmonics become more 
significant. 



I. INTRODUCTION 

Ground-based gravitational wave detectors like LIGO 
and Virgo will soon detect the complicated gravitational- 
wave signature of merging binary compact objects. In 
general, all aspects of that signal will be significantly 
modulated by each compact binary's spin: misaligned 
spins break symmetry in the orbital plane, generally lead- 
ing to strong modulations in the inspiral pQ, merger, 
and ringdown. In the early stages, these modulations 
are well-approximated by quasistationary radiation from 
a circular orbit, slowly rotating and modulated as the 
orbital plane precesses [DBJ- At late times, numerical 
relativity simulations of mergers can efficiently predict 
gravitational radiation from a variety of spins and or- 
bits. A well-chosen static or time-evolving frame makes 
it easier to compare their outputs to one another and to 
these analytic models. Several methods to choose a frame 
have been proposed, using two broadly distinct algo- 
rithms. One proposed method chooses a frame such that 
the instantaneous (l,m) — (2,2) emission is maximized 
[H [5] . Alternatively, a preferred frame follows by averag- 
ing a tensor constructed from the action of the rotation 
group on asymptotic emission over all orientations [BJ. 
Though originally investigated numerically, both meth- 
ods are analytically tractable and can be applied to ex- 
isting tables of post-Newtonian waveforms' spin-weighted 
spherical harmonic modes [3J. In this paper we evaluate 
(£( a £f,)) in the post-Newtonian limit. 



II. PREFERRED ORIENTATIONS AND 
QUASICIRCULAR ORBITS 

We explore gravitational wave emission from adiabatic 
quasicircular binary inspirals. These orbits can be pa- 
rameterized by the two component masses mi, m 2 ; com- 
ponent spins Si, S2; orbital separation vector r; and re- 



duced orbital velocity v = d t r. Following the notation 
of [3J, we perform our post-Newtonian expansions using 
dimensionless mass and spin variables: 
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S = x-+ Sx+ 

M \mi m 2 

We note that S is essentially a dimensionless version 
of the spin variable A from [7], as the are related by 
— M 2 S = A. To describe the orbit and its orien- 
tation, we define the "Newtonian" angular momentum 
Lyv = Mrjr x r; the unit separation vector h = r/|r|; 
and the unit velocity vector v — v/|v|. 

The tensor (Ci a Cb)j is defined by averaging a certain 
angular derivatives of the Weyl scalar ^4 over all ori- 
entations, then normalizing the result [6J. Because this 
tensor is constructed from rotation group generators C a , 
this normalized orientation average can be calculated ef- 
ficiently from spin-weighted spherical harmonic basis co- 
efficients ipii m - 



(£(aA)) = 



/ dntp 4 * (t)c {a c b} tp 4 (t) 



j<m\vb 4 \ 2 

Ew j^m^hlm ( lm ' \C (a C b) \lm) 
J dn\yj 4 \ 2 



(5) 



'Electronic address: evano@gravity.phys.uwm.edu 
^Electronic address: oshaughn@gravity.phys.uwm.edu 



where in the second line we expand vb 4 = 
T,i m ' t l ; 4i m {t)Y l '^ 2 \6,<j>) and perform the angular 
integral. In this way, the average tensor ^£( a £ b )) at 
each time t can be calculated algebraically; see the 
appendix for explicit formulae. To better understand 
this average, we will contrast this expression with a 
similar average, |£( a £ & )}, defined by replacing ip 4 — > h 
in the above expression. 

Several general conclusions about the preferred ori- 
entations implied by (£( a £^)) follow from the defini- 
tion. First and foremost, as gravitational radiation is 



2 



quadrupolar and predominantly driven by mass moments 
during inspiral, the tensor should be well-approximated 
by the corresponding limit implied by instantaneous, 
purely quadrupolar emission [5J: 



(£(aA)) 



LNa^Nb {m 2 ) 

(l(l + l)-m 2 ) 



[h a n b + v a v b ) (6) 



To the extent that the waveform is dominated by the 
leading-order (I, \m\) = (2, 2) modes, this tensor is diago- 
nal, with eigenvalues (4, 1, 1). The dominant eigendircc- 
tion V of (£( a £{,)) is therefore nearly the "Newtonian" 
orbital angular momentum. Second, from perturbation 
theory, the dominant eigendirection V can be modified 
at leading order only through terms of the form L a X b for 
some vector X. 

A post-Newtonian expansion of (Cr a C b \) has addi- 
tional convenient structure when (£(- £(,0 is restricted 
to a constant-/ subspace. Specifically, in Eq. ([5]), we 
employ terms in the numerator and denominator from a 
single I. Because only terms of constant I enter, such an 
average has fixed trace: 



Tr(£ (a £ b) ) / = 1(1 + 1) 



(7) 



A series expansion of (£(„£&)) therefore consists of an 

expected zeroth-order term [7"'"], followed by symmetric 
tracefree tensors [yl^*^] of increasing order and tightly 
constrained symmetry: 



= t« 
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k=0 



= [l I L Na L Nb + l/2{h a n b + v a v b )\ (9) 

We seek to understand what tensors arise in a pertur- 
bation expansion and to determine how they impact the 
preferred orientation V . 



III. ROTATION TENSOR IN THE 
POST-NEWTONIAN LIMIT 



We employ a post-Newtonian model for the spin- 
weighted spherical harmonic modes of either the Weyl 
scalar, tp4, or the strain, h, for a generic precessing bi- 
nary to construct the tensor (Cr a C b \) or |£( a £ b )}. The 
modes hi m are provided in [3] in a specific coordinate 
frame, where the z-axis is aligned with Jo, the total an- 
gular momentum at some initial time. In this frame, the 
instantaneous "Newtonian" orbital angular momentum 
Ln = Mrjr X f has the form 



zcos(l) + sint(cosai + sin ay) (10) 



We find that our expressions for the rotation tensor are 
simplest when given in a frame aligned with the Newto- 



nian orbital angular momentum in which 

n = {cos $ , sin $ , 0} 
v = {-sin<I>, cos$, 0} 
L N = {0,0,1}. 
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The modes of [3] will be given in this frame simply by 
evaluating the general expresions for i — 0, a — it. We 
derive our expressions for (Lr a C b \\ and {£( a £;,)} in this 
frame, but then express them in a frame-independent 
matter in terms of the vectors n, v and L^. Thus, 
(£( a £f,)) and {£( a £{,)} can be given in any frame simply 

by finding the components of h, v and in that frame. 

In the non-precessing case, we find that the PN cor- 
rections to (£( a £ & )) and {£( a £f,)} (and their restrictions 
to constant I subspaces) can always be expressed as lin- 
ear combinations of three symmetric, trace-free matri- 
ces. These matrices (and their components in the frame 
aligned with Ljv) are: 
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A. Weyl scalar (il>i) in the Post-Newtonian 
quasistationary limit 

To evaluate (£( a £;,)), we convert from the radiated 
strain moments hi m provided in [3J to ipi m by taking time 
derivatives: 
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and organizing this expression self-consistently in powers 
of v. 

For non-precessing binaries, all derivatives save $ and 
v are zero. At leading order, $ = w or b = v 3 /M and 
v = (32r]/5M)v 9 , with corrections to each of these 
known to 3.5PN relative order. Also note that the spin- 
independent terms in hi m are known to v 6 (3PN) [5J H] 
beyond leading order and the spin-dependent terms to v A 
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order (2 PN) [3). Because v is order v 6 higher than 4>, v 
can be neglected so long as one is interested in expressions 
for tpi m and (£(„£&)) accurate to less than 3PN order. In 
particular, we need to consider v only for 3PN-accurate 
ipim for non-spinning binaries. 

For precessing binaries, the hi m are known only to 
1.5PN order [3J, so this is the highest order to which we 
can obtain ipi m and (Cr a L^\). Spin effects enter hi m at 
1PN and 1.5PN order and $ at 1.5PN order. We must 
also account for the derivatives of hi m with respect to 
the spins, because Si oc v 5 which is 1PN order relative to 
<I> HH] ■ Lastly, we must also keep track of I and a (or 

equivalently Ln, as the angles i and a describe the orien- 
tation of Ln), since Ljv = Ql n x Ln oc v 6 (see Eq. (43) 
below) is 1.5PN order relative to $. 

For geometrical reasons, the corrections from the 
derivatives a and I depend only on £Il n , the vector 
around which Ljv precesses, specifically as the combi- 
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Note that the motion of Ljy creates a 1.5PN correction 
(relative to $) in ^»j m , while we have defined T to be 
dimensionless by canceling out the v 6 /M dependence of 



B. Rotation tensor for non-precessing binaries 

For a non-precessing binary, the tensors (£( a £f>)) and 
{£(„£{,)} have the following frame-independent expres- 
sions: 
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Because gravitational wave emission is predominantly 
quadrupolar, these tensors are T^> at leading order, with 
PN corrections appearing at 1PN order. 

We can also consider these tensors restricted to sub- 
sets of modes with constant Z, consistently in both the 



numerator and denominator of the average. To 1.5PN 
order, when we restrict these tensors to the I = 2, 3, 4 
subspaces we find 
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Note that in all cases these tensors are block diagonal, 
with a block corresponding to the orbital plane spanned 
by h and v and a block for L^- In particular, there are 
no off-diagonal terms to couple the contributions along 
Ljy to those in the orbital plane. This means that the 
dominant eigenvector will remain along Ljv at all times 
with an eigenvalue ~ I 2 and there will be two eigenvec- 
tors in the orbital plane with eigenvalues 1/2. While the 
eigenvectors will be constant, the PN corrections will in- 
troduce small, time-dependent variations in the eigenval- 
ues. In fact, for nonprecessing binaries this direction is 
tightly enforced by symmetry and cannot be modified, 
for example, by omitted nonlinear memory terms; see 
Appendix [A"| for details. 

The tensors extracted from constant-^ subspaces can 
behave pathologically when multipoles of that order are 
strongly suppressed by symmetry. In particular, the 
(3, 3) and (3, 1) modes vanish in the equal mass case 
where S — 0. As a result, the expressions for (£( a £f,)) 3 
and {£( a £f>)} 3 will diverge in the limit S — > because 
they contain a factor (1 — 4ry) = S 2 in the denominator. 



are no longer block-diagonal. These off-diagonal terms 
depend on the transverse components of spin and mix 
contributions in the orbital plane with those along Lm- 
As a result, the eigendirections will vary in time as the 
binary precesses and they will not lie along Ljy and in 
orbital plane. 



As in the non-precessing case, we derive the expressions 
for ^£( a £ b )) and {£( £;,)} in a frame aligned with the in- 
stantaneous L n but express them in a frame- independent 
way so that they can be used in any frame. As in the non- 
precessing case, all of the tensors can be expressed as lin- 
ear combinations of a small collection of matrices. The 
new matrices needed for the precessing case (and their 
components in the frame aligned with the instantaneous 
Ljv) are: 



C. Rotation tensor for precessing binaries 

In the precessing case, when the spins are not aligned 
with the orbital angular momentum, the rotation tensors 
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and also Mac 1 X+cLno and M^} X+cLNb, which arc the 
equivalent to the expressions above, but with the com- 
ponents of x+ replacing those of S. 



The and tensors and their constant 

I restrictions for precessing binaries are given by 



5 



(c {a c b) y 



{£(<»£&)} 

( C (a C b)) 



I „ „ \ NP a 

(£ (o A)) +w 3 
+ 6u 3 L Ar(a T b) 

Ptcc r - „ I NP 3 

prcc / - \ NP 3 

+ 6w 3 L A r (a T b) 

NP 



{ £ («A)}3 
(£(aA)} 



prcc 
4 



(£( a A)) 4 



581 



M« — <5S + 29r/ X+ L 



96 



41 



6 -55 + 29t ? x+J L Nb) -Mgl i^fSS, + 21,,\_,. j I A ,„ 
581 



35 



M&>(— JS + 29»7X+) ^6) + -M ( ^ < S c i JV6) 



96 (° 



M^ 1 



np ,7680(1- 3t?) (1) 



(r ( „A)) 3 = (AaA)) 3 6151(1 _ 4?? ) M (^+^) 

= {£ {a £ b) } 3 " v 19(1 _ 4t?) M{ a > cX +cL Nb) 



- ( £ (aA)) 



NP 

4 

NP 



(31) 
(32) 

(33) 
(34) 

(35) 

(36) 

(37) 
(38) 



At high multipole order, the accuracy of these ex- 
pressions is inevitably limited by the accuracy of post- 
Newtonian expansions, which contain terms explicitly 
coupling to spins only to modes I < 3 [10]. Hence, 
the I — 4 subspace does not have any off-diagonal terms 
through 1.5PN order, though it will presumably acquire 
such terms at a higher PN order. 

By examining Eq. (4.17) of [3], we can understand the 
origin of these off-diagonal terms with transverse spins. 
It may be surprising that they do not appear until 1.5PN 
order, since h 2 2 has such terms already at 1PN order. 
However, as explained in Appendix [AJ these transverse 
terms are the result of coupling modes with m values 
that differ by 1. As a result, the 1PN transverse spin 
terms in h 22 must couple with the leading 0.5PN order 
term in h 2 \, so that the off block-diagonal contribution 
to {£( a £ b )} is 1.5PN order. These terms in h 22 and h 2 i 
are both proportional to the first harmonic of the orbital 
phase, so this orbital scale dependence cancels out and 
creates a contribution to {£( £ b ) } of the form of Eq. ( 29 ) 



(i.e. M£>S c L Nb for some combination of spins S) . There 
are also 1.5PN order transverse spin terms in h 2 \ propor- 
tional to the second and zeroth harmonic that couple to 
the leading order term in h 22 . These create contribu- 



tions to {£( a £ b )} of the form of Eqs. (29) (M$) S c L Nb ) 



and (30 1 (M ac S c LNb)> respectively. 

In a similar manner, /i 33 and /131 have transverse spin 
terms proportional to the second harmonic which couple 
to the leading term in /132 an d cr eate off-diagonal terms in 
|£( Q £ b ) } of the form of Eq. S. These terms are 2.5PN 
order, and so they do not appear in Eq. (32). However, 
because the leading I = 3 terms are 0.5PN order, and 



in computing j£( a £ b )} 3 we restrict both the numerator 
and the denominator to I = 3, these terms are 1.5PN 
relative to the leading order of {£( a £ b )} 3 and therefore 
appear at that order in Eq. ( [36] ) . The I — 4 modes do 
not have transverse spin terms through 1.5PN order, and 
so we do not see any off-diagonal terms in {£( a £ b )} 4 . 



D. Broken symmetry from precession 

For a precessing binary at 1.5PN order, no 
longer has the Newtonian, orbital, or even total angu- 
lar momentum as an eigenvector. The orbital and total 
angular momenta are given by the spin-dependent ex- 
pressions [331 02]: 
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where we have explicitly expanded L to 1.5PN order. 
Because the Newtonian angular momentum at leading 
order is cx u -1 , the term L50 m this expansion is v 3 
beyond leading order. 

At and below 1PN order L = L N . The 1.5PN or- 
der spin-orbit contribution to the orbital angular mo- 
mentum, L50 breaks this alignment and introduces a 
component in the instantaneous orbital plane that oscil- 
lates at twice the orbital frequency. Likewise, the spin- 
dependent terms at order v 3 in (£( a £f,)) and {£( a £b)} 
are not block diagonal and time dependent. As we will 
show by concrete example in Section |IV[ these two time- 
dependent expressions do not conspire to evolve together: 
L is not an eigenvector of either (£( a £;,))j {£(a£b)}, or 
most (but not all) of the tensors produced from constant-/ 
subspaces. To demonstrate this explicitly, one can labo- 
riously evaluate the expressions 
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An eigenvector X a of the three-dimensional matrix A a b 
necessarily has t a bcX h ' A c( iX d — \e a b c X b X c = 0. For 
(C( a Cb)) and {£( a £b)}, we find each of these vectors are 
nonzero atl.5PN order, once the off-diagonal precession 
terms enter. 



E. Memory terms 

In the calculations above, we use the "ready to use" 
waveforms provided by Arun et al. [3] . These expressions 
do not include nonlinear memory terms, which enter at 
leading PN order and beyond. In terms of the modes 
used in this paper, these terms would produce a strain 
multipole h 2 o comparable in magnitude to h 22 for non- 
precessing binaries [131 H] • 

No complete expression exists in the literature for non- 
linear memory for generic spinning systems. Memory 
terms depend on the integrated past history of the bi- 
nary; they therefore can differ in magnitude, angular de- 
pendence, and functional form versus time depending on 
how the spins have evolved in the past. That said, we 
anticipate our estimate of this particular preferred ori- 
entation (i.e., the dominant eigenvector of (C( a Cf,)}) will 
not change significantly when memory is included. For 
example, because memory terms cannot modify the re- 
flection symmetry for nonprecessing binaries, they like- 
wise cannot change the preferred orientation away from 
the normal to the orbital plane. 



IV. NUMERICAL EVALUATION AND 
DISCUSSION OF RESULTS 



{£( a £f,)}, differ from L, j, and and Ln, with a differ- 
ence that grows with time. To illustrate this difference, 
we have simulated the evolution of two precessing, inspi- 
ralling binaries: a 10 -I- 1.4M BH-NS binary where the 
black hole has a maximal spin tilted 30° relative to the 
Newtonian orbital angular momentum at its initial fre- 
quency of 40 Hz and a 60 + 30A/ Q BH-BH binary where 
both black hole spins are maximal and in the orbital 
plane with a 90° angle between the two spins at an initial 
frequency of 10 Hz. Figs. [I] and [2] plot the trajectories of 
Ljf, L, J, V and Vh for these two binaries. 

The equations of motion were numerically integrated 
using the lalsimulation 1 implementation of the adi- 
abatic, spinning Taylor T4 approximation described in 
[16] . Specifically, this code evolves the binary phase and 
frequency, spin and Newtonian orbital angular momen- 
tum vectors equations according to [see, e.g., [II IIB"! 117]]: 
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where E and T are the PN-expanded binary energy 
and gravitational-wave flux (with spin contributions in- 
cluded), respectively. F /(MdE/dv) is then re-expanded 
as a polynomial in v, as described in [18] . At each time, 
we obtain values for (£( a £b)) and \C( a Ct,) ) by taking the 
dynamical variables computed via Eqs. (43 ) and plugging 



these into Eqs. (31 )-(38 1. We then compute the principal 
eigenvectors and eigenvalues of (£( a £b)) and {£( a £b)} 
(and their constant-/ subspaces), which are used to cre- 
ate the plots in this section. 



Our calculations above have demonstrated that V 
and Vh, the principal eigendirections of an d 



1 The lalsimulation package encapsulates and standardizes sev- 
eral waveform generation codes; see the LIGO Algorithms Li- 
brary 1151 for more details. 




FIG. 1: Comparing trajectories 1: BH-NS binary: 

These panels show the trajectories of Ljv (black), L (green), 
J (cyan), V the principal eigendirection of (C( a C b )) (blue) 

and Vh the principal eigendirection of {£( a £j,)} (red) for a 
10 + IAMq BH-NS binary where the black hole has maxi- 
mal spin tilted 30° relative to the Newtonian orbital angular 
momentum at an initial frequency of 40 Hz. The frame is 
chosen so that the z-axis coincides with the initial value of 
J. The top panel shows the first precessional cycle, which 
spans a frequency range of roughly 40 — 47 Hz. The bot- 
tom panel shows a very late precessional cycle shortly before 
merger would occur, covering a frequency band of 200 — 600 
Hz. 



A. Orientations 

The trajectories shown in Figures [l] and [5] demonstrate 
that, for these two prototype binaries, the four direc- 
tions L,Ln,V and Vh generally do not coincide. Though 
these four orientations have very roughly the same fea- 
tures - all precess - all four exhibit distinctly different 
trends. For example, the Newtonian angular momentum 
Ln (black) precesses smoothly for our analytic model. 
The orbital angular momentum L oscillates (nutates) 
around Ljv- The preferred orientation V extracted from 
the Weyl scalar (blue) also precesses smoothly in our co- 
ordinates, tracking with but offset somewhat from Ljy. 
Finally, the preferred orientation Vh extracted from the 



FIG. 2: Comparing trajectories 2: BH-BH binary: This 
plot is the same as Fig. [I] but for a 6O + 3OM0 BH-BH binary 
with both spins maximal, perpendicular to the Newtonian 
orbital angular momentum and perpendicular to each other 
at an initial frequency of 10 Hz. This plot shows the entire 
inspiral of the binary, reaching a GW frequency of about 85 
Hz at the end. 



strain (red) oscillates around V. 

As expected from their functional form - a power series 
in v - the differences between these four quantities grow 
more pronounced as v grows, closer to the merger event. 
As a consequence, in the fixed sensitive band of ground- 
based detectors, higher-mass sources exhibit more pro- 
nounced differences between these four orientations. We 
cannot reliably extrapolate to the late stages of merger. 
However, as suggested by previous studies [15], we expect 
substantial differences can accumulate between these ori- 
entations. 

Among these four orientations, only and V evolve 
smoothly, without exhibiting strong oscillations on the 
orbital period. By contrast, the direction Vh extracted 
from the strain oscillates significantly on the orbital pe- 
riod. The differing behavior of Vh and V can be under- 
stood entirely from the magnitudes of the terms of the 
form of Eq. (29) and Eq. (301 which appear in Eqs. (31) 



and ( 32 ) . Note that the former depend only on the spin 



components, and so create smooth perturbations to the 
eigenvectors which vary on the precessional time scale. 
The latter vary as twice the orbital phase, and so create 
perturbations that oscillate much more rapidly. 



So, by comparing Eqs. (31) and ( |32[ ) we see that the 
coefficients of M^' in (£( a £ b )) and |£( a £ b )} have very 

similar magnitudes, and so offset V and Vh from Ljv by 
similar amounts. However, the coefficient multiplying 
in (CuCfy) is always much smaller than the co- 
efficient multiplying M^t '. Therefore, these oscillations 
are strongly suppressed and are not apparent in Figures 
[T] and [2] On the other hand, the coefficient multiply- 




Time (s) 



03 -0.2 
Time (s) 



FIG. 3: Comparing eigenvalues: The fractional deviation (A/A — 1) in the dominant eigenvalues (where A is the unperturbed 
eigenvalue and A is the actual eigenvalue) of (top panels) and {£( a £f,)} ; (bottom panels), for the BH-NS binary orbit 

described in Figure[T] The blue curve indicates I = 2; the red dashed curve shows I = 3. For reference, a black solid line is shown 
at 0. Left panels: Eigenvalues of (jCi a Cty} (top) and {£( a £{,)} (bottom) for the full evolution shown. Right panels: As above, 
but limited to a short interval immediately preceding the last simulated post-Newtonian orbit, where the minimum-energy 
condition breaks down. 



ing M^ b ' in |£( a £ b ) } is nearly as large as the coefficient 

multiplying . Therefore, the orbital-scale oscilla- 
tions are much more significant and creates the distinct 
"crown" shape in the trajectory of Vh in Figures [l] and 



vanish if and only if S — 0. But this can only happen if 
Si /mi = S2/7712 which would mean that the dimension- 
less spin of the larger body must be extremely tiny and 
the system will again not precess significantly. In inter- 
mediate cases, the off-diagonal terms would vanish only if 
the contributions from the x+ symmetric spin terms and 
the S antisymmetric spin terms were carefully tuned to 
cancel out. However, one cannot simultaneously ca nce l 
derivatives oi L N , where T is given by Eq. (|20j). One can the terms multiplying both M$ and in Eqs. (E) 



Another difference between Vh and V is the presence 



of the L N r a T b \ term in Eq. (31 1 that arises from the time 



see that T will always point back towards the center of 
the precessional cone (J), so that this term always tends 
to shrink the precessional cone of V relative to Vh- In 
particular, in Figs. [T]and[2] we see that the trajectory of V 
lies just above the oscillations of the trajectory of Vh and 
we have checked that if the L^^T ^ term were removed 

V would pass through the center of the oscillations of Vh- 

Let us consider V and Vh in the equal mass limit 
(5 -> 0, r) -> 1/4). From Eqs. pi]) and Q it is clear 



and (32). Therefore, we can conclude that for any pre- 



the off-diagonal terms will vanish if and only if x+ = 0, 
which means the spins must be equal and opposite. But 
in this case the total spin vanishes and the binary does 
not precess. On the other hand, in the extreme mass 
ratio limit (<$—»• 1, 77 0) the off-diagonal terms will 



cessing binary the dominant eigendirections V and Vh 
cannot coincide with Ljq. 

To this point we have described the preferred orienta- 
tions V and Vh associated with (C/ a Cf,)) and {£( a £b)}, 
summing over the contribution from all subspaces. We 
can also define the two preferred orientations associated 
with the constant-/ subspaces of and {£( a £;,)}. 

One would expect (and Figure |i] confirms) that, because 
leading-order quadrupole emission dominates the gravi- 
tational wave signal, V2 {Via) would be extremely close 
to V (Vk)- By contrast, Figure [4] also demonstrates 
that Vh3 — Ln. To the PN order considered here, 
V4 = V M = L N , since (£( Q £ h )) 4 and {£( Q £ b) } 4 do not 
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corresponding to the orbital plane and the Newtonian 
orbital angular momentum, such as those in Eqs. (21) - 



(28), can either raise or lower the eigenvalue depending 



-2a -1.5 
Time {s) 



FIG. 



4: Orientations of subspace eigenvectors: Top 

A plot of 1 - Vi ■ V (solid blue) and 1 - V h2 ■ V h (dashed 
red) versus time for the BH-BH binary orbit described in Fig- 
ure [2] Bottom panel: A plot of 1 — V3 • Ljv (solid blue) and 
1 — Vh3 ■ Ln (dashed red) versus time for the same BH-BH 
binary. The eigenvectors from the I = 2 (I = 3) subspaces 
are rather close to the eigenvector of the full rotation tensor 
(L N ). 



have off-diagonal terms. 



B. Eigenvalues 

Generally, the eigenvalues of (£( £{,)) are powerful di- 
agnostics for gravitational wave signals with multiple har- 
monics, particularly in the presence of precession. For ex- 
ample, to an excellent approximation, the eigenvalues of 
any constant-/ subspace orientation tensor (£( a £ & )) ; are 
I 2 /2 and 1/2. As noted in [5J, the dominant eigenvalue 
of (£( a £ft)), summing over many constant-/ subspaces, 
scales as (m 2 ^ and therefore reflects the relative impact 
of higher-order harmonics to the signal. 

To demonstrate the value of these eigenvalues as a di- 
agnostic, in Figure [3] we show the relative difference in 
the eigenvalues of (£( a £ft)) and {C( a Cb)} from the naive 
leading-order predictions described above. The PN cor- 
rections create rather small perturbations to the eigen- 
values. The PN corrections within the diagonal blocks 



on their sign. The off-diagonal PN corrections will always 
increase the dominant eigenvalue, regardless of their sign. 
See Appendix |B] for details. 

The PN corrections of course become largest at late 
times, when v increases rapidly, and this is why devia- 
tions in the eigenvalues grow rapidly at the end. For the 
/ = 3 subspace, the PN correction in the Ljv block [i.e. 
multiplying in Eqs. (25) and (26)] is always neg- 



ative and hence the dominant eigenvalue is always less 
than its leading-order value. For the / = 2 mode, note 
that the 1PN and 1.5PN corrections have the opposite 
sign. At early times, the 1PN term dominates and tends 
to decrease the dominant eigenvalue. Depending on the 
values of 5 and S, the 1.5PN correction can actually be- 
come larger than the 1PN correction and change the sign 
of the perturbation to the dominant eigenvalue (as hap- 
pens in our BH-NS binary). 

As with the trajectory V, the eigenvalue of (£( a £ b )) 2 
and ^£( Q £ h )) 3 changes smoothly on the precession and 

inspiral timescales. By contrast, like the trajectory Vh, 
the eigenvalues of {£( a £ b )} 2 oscillate rapidly on the or- 
bital timescale; see the bottom two panels. [At the post- 
Newtonian orders available at present, the orientations 
and eigenvalues of {C^ a C b - ) } 3 do not oscillate signifi- 
cantly.] 



V. NUMERICAL EVALUATION AND 
SUBTLETIES OF SELF-CONSISTENT 
POST-NEWTONIAN EXPANSIONS 

The post-Newtonian series converges slowly in powers 
of v , itself only between 0.1 — 0.5 for our examples. The 
details of how PN expressions are truncated or resummed 
can quantitatively (though not qualitatively) change the 
conclusions described above. Our concern is not aca- 
demic: there exist several straightforward alternatives 
to this calculation which will cause (£( a £{,)), {£( a £ b )} 
and the precession equations to effectively differ by un- 
known, higher PN order error terms. In this section, 
we describe how these subtle issues associated with the 
post-Newtonian series can impact our results. 



Numerical evaluation of preferred orientation 
from mode coefficients 



To compute (£( a £{,)) and {£( a £f,)}, we consistently 
expand both the numerator and denominator in a post- 
Newtonian series. Such a rational function form would 
be a completely reasonable choice to represent the rota- 
tion tensors. Instead, for simplicity in the expressions 
and their interpretation, we have chosen to re-expand 
the rational function as a polynomial via a Taylor series. 
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The two agree to 1.5PN order, but differ by some un- 
known term at 2PN order. This is exactly analogous to 
the difference in TaylorTl and TaylorT4 waveforms de- 
scribed in |18j . with the rational form corresponding to 
TaylorTl and our polynomial expressions corresponding 
to TaylorT4. 

Furthermore, a straightforward numerical implemen- 
tation to generate {£( a £{,) } will also differ slightly from 
the rational and polynomial expressions described above. 
This method would employ Arun et al. [3] to compute 
%i m , then apply the formulae in Appendix [A] to com- 
pute {£( a £;,)} directly. However, as {£( a £{,)j involves 
ratios of products oc h 2 , this direct calculation implicitly 
includes many higher-order post-Newtonian terms. For 
example, in these products a v 3 correction to one hi m 
could multiply a v 3 correction to another hi m i to give a 
contribution to (£( £(,)) oc v 6 (3PN). This is only a par- 
tial PN correction, of course, because if the hi m 's were 
complete to 3PN there would be many more terms at v e . 

Because of these subtleties, the values of (£( a £f,)) 
and {£( a £;,) } computed via polynomial, rational func- 
tion and direct numerical computation will differ from 
one another. While all of the plots in this paper were 
computed using the polynomial expressions of Eqs. (31 1- 



( 38 ), we have checked that the rational and direct numer- 



ical computations are quite similar. There are quantita- 
tive differences, but not qualitative ones. The differences 
among V, Vh, L and Ln are quite similar, the trajec- 
tory of Vh still has its distinctive "crown" shape, and our 
conclusions are not changed. 



B. Orbit-averaging the precession of Ln and L 

Note that to create the plots in Figures [I] and [2] we have 
used an orbi t-ave raged formula to describe the precession 



of L N , Eq. (43c). As a result, L N moves on a smooth 
precession cone and L nutates on an orbital timescale 
because of its Lso term. This is in sharp contrast with 
Fig. 6 of j2], who find that L varies smoothly, while 
Ln nutates on an orbital timescale. This is presumably 
because the authors of [5] used an orbit-averaged preces- 
sion equation for L (rather than Ln). As explained in 
the text immediately after Eq. (4.17) of [7], Ln and 
L have the same orbit-averaged precession equation, be- 
cause Ti so — when averaged over a circular orbit. 

We emphasize that our expressions for (^C^ a C^ and 
|£( a £ft)} in Sec. Ill do not make any sort of orbit- 
averaged approximation. Indeed, the use of a PN expan- 
sion for hi m and ipiim * s the only approximation involved 
in deriving these expressions. The orbit-averaging ap- 
proximation is applied when integrating the equations of 
motion to generate time series for Ljv and the spin vec- 
tors to plug into the expressions. Thus, the details of the 
plots in Sec. |IV| may vary depending on if or how orbit 
averaging is applied. 

Therefore, when constructing PN waveforms or per- 



forming further studies to find preferred directions for 
precessing binaries, one should carefully consider how the 
precession equations handle this orbit averaging. Several 
authors have investigated more accurate formulations of 
the precession equations which do not assume orbit av- 
eraging [7j [17l l20Tf2"2] . Given the stringent accuracy re- 
quirements needed for parameter estimation, comparison 
with numerical simulations and other applications, we ex- 
pect that the difference between non-orbit-averaged and 
orbit-averaged precession equations could be significant. 
Therefore, we recommend future studies incorporate ana- 
lytic waveforms with non-orbit-averaged precession equa- 
tions to be as accurate as possible. 



VI. CONCLUSIONS 

Gravitational wave emission from precessing compact 
bianries is expected to appear simpler in a suitably- 
chosen frame. To leading post-Newtonian order, this 
frame should correspond to the orbital angular momen- 
tum at large separations. Late in the inspiral, many 
plausible reference orientations exist. In this paper, we 
have demonstrated that one choice for that orientation 
(the prinicipal eigendirection of (£( a £{,))) does not agree 
with cither of two others (the orbital and total angular 
momenta) . 

The suitability of this or any other reference frame is 
best judged by whether it facilitates useful insights and 
calculations. In this paper we demonstrate the tensor 
(£( a £f,)) provides a particularly tractable way to extract 
preferred orientations from asymptotic radiation. Our 
approach allows different preferred orientations to exist 
at each multipolar order. Previous papers have demon- 
strated this method useful in analyzing numerical relativ- 
ity simulations [SJ [05] Finally, in the text and appendix 
we provide methods to evaluate this average using stan- 
dard tables of hi m modes. 

Our investigations suggest the proposed orientations 
will differ from one another for any precessing binary, 
regardless of its mass ratio or spin orientation. Further- 
more, we note that in the PN regime the behavior of 
these preferred directions can be noticeably affected by 
the details of how the PN series is truncated or resummed 
and also by the manner in which orbit averaging of the 
precession equations is performed. This suggests that 
precessing PN models could benefit greatly from using 
more general precession equations which do not assume 
orbit averaging. 
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Appendix A: Evaluating the average 

For reference, we provide explicit formula that relate 
the mode amplitudes ipi m of the Weyl scalar to the 
orientation-averaged expression (£( a £ b )). We first cal- 
culate two real and two complex quantities: 



h = 



h 



- (ip,L + L + ip) 



W>,L + 0L 2 + 1/2)V) 

Y Ci m (m + l/2)lpl m+1 l/>i r , 



I 111 



= ^[/a + l)-m 2 ]|^ m | 2 

lm 

= (^,L z L z lfj) = y^TO 2 |^; m | 



(Al) 
(A2) 

(A3) 
(A4) 



where c; m = + 1) — m(m + 1). In terms of these 
expressions, the orientation-averaged tensor is 



J + Re(J 2 ) lml 2 ReJi 
/ -Rc(/ 2 ) ImJx 

I zz 



(A5) 



turn (oc z) will be the principal axis of (£(„£;,)) if and 
only if ii = 0. The prefactor C/ m = c/ m (m + 1/2) in 
the definition of I\ is antisymmetric about m = — 1/2 - 
that is, under a transformation to — > — (m + 1). There- 
fore, one way to make 1\ — and thus L a prin- 
cipal axis of (£( a £f,)) occurs if ">P* m +i'4 > im is symmet- 
ric about to = —1/2. For nonprecessing binaries the 



strain modes satisfy hi 



(reflection symmetry 



In a frame aligned with L, the orbital angular momcn- 



through the xy plane) and evolve in phase with the or- 
bit (him = e 4m *^4;,|m| for some real Ai m ). Therefore, for 
each constant-^ subspace, term-by-term cancellation oc- 
curs in ii , independent of the details of A\ m As a concrete 
example, we expand out I\ for I = 2 as 

h = e^[A 20 A 21 (C2-i+C 2 o) + A 21 A 22 (C 2 -2 + C 2A )} 
= (A6) 



Appendix B: Perturbations of the eigensystems 

First, let us consider a non-precessing binary and study 
(£( Q £b)) ; in a frame aligned with Ljy. From Eq. j9|), wc 
see that to leading-order 

(1/2 0\ 
(£(.A)>, - T {1) = 1/2 , (Bl) 
\ 0/7 

and of course the eigenvalues are I 2 , 1/2 and Z/2. The 
PN corrections are small perturbations to the tensor, so 
when we add them in it takes the form 



(C {a C b) ) l = T W + eiAf« + e 2 M( 2 ) + e 3 M^ 



'1/2 + a + e 2 cos 2$ 

e 2 sin2$ 1/2 




e 2 sin 2$ 
- e\ — e 2 cos 2$ 
l 2 + e 3j 



(B2) 



In this case, the eigenvalues become I 2 + £3 and 1/2 + 
ei±e 2 . So, PN corrections to the diagonal blocks (corre- 
sponding to either the orbital plane or the Newtonian or- 
bital angular momentum) will linearly perturb the eigen- 
values within that block. They do not change the dom- 
inant eigenvector. While they do weakly break degener- 
acy in the orbital plane, setting two distinct eigenvectors 
rotating in the orbital plane, they do so weakly; the space 
spanned by these two vectors is still the subspace of the 
unperturbed system. 

Now let us consider the perturbations from the off- 
diagonal terms that are present in precessing binaries. 
For simplicity, we will consider the / = 2 case, although 
the same argument can be applied to any value of I and 
it is straightforward to derive more general expressions. 
In this case, in a frame aligned with the instantaneous 



ZrjVj our perturbed rotation tensor will have the form 

/l eA 

{C (a C b )) 2 = 1 e a . (B3) 
\ei £2 4/ 

The eigenvalues of this system are 

1(5 + 3^1 + ^ + 61)) ~ (B4) 

^(5 - ZsJl + \(e\ + 4)) * 1 + ! ^ (B5) 
1 (B6) 

In particular, we see that the presence of these off- 
diagonal terms will always increase the magnitude of the 
dominant eigenvalue, regardless of the sign of these per- 
turbations. It is also worth noting that the dominant 
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(unnormalized) eigenvector is given by 



V 2 = L,e 2 l ll + Jl + ^4 + 4) 



(B7) 



fact that oscillatory off-diagonal terms such as those of 
the form in Eq. ( 30 ) are directly responsible for oscilla- 



tions in the principal eigendirections. 



so that V2x (^2j/) is linearly proportional to (C,( x C z y} 2 
((£(j / £ z )) 2 ). This shows quite clearly, for example, the 
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